Microscopic Theory of Rashba Interaction in Magnetic Metal 
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Theory of Rashba spin-orbit coupling in magnetic metals is worked out from microscopic Hamilto- 
nian describing d-orbitals. When structural inversion symmetry is broken, electron hopping between 
d-orbitals generates chiral ordering of orbital angular momentum, which combines with atomic spin- 
orbit coupling to result in the Rashba interaction. Rashba parameter characterizing the interaction 
is band-specific, even reversing its sign from band to band. Large enhancement of the Rashba pa- 
rameter found in recent experiments is attributed to the orbital mixing of 3d magnetic atoms with 
non-magnetic heavy elements as we demonstrate by first-principles and tight-binding calculations. 

PACS numbers: 75.60.Jk 



Control of the local magnetization direction in a metal- 
lic ferromagnet through delivery of torque by the current- 
carrying electrons is one of the major endeavors of spin- 
tronics community nowadays [1]. Prompted by the the- 
oretical discovery of a new type of spin transfer torque 
(STT) arising in Rashba-coupled bands [2-5] and subse- 
quent experimental indication thereof [6] , several propos- 
als [7, 8] were made recently to uncover a complete set 
of STT's permissible in ferromagncts with the Rashba- 
coupled bands, or Rashba ferromagncts for short. Exper- 
iments to measure Rashba-induced STT were performed 
for a very thin atomic layer of magnetic atoms with sev- 
eral layers of heavy atoms such as Pt grown on top [6]. 
It was predicted that the Rashba interaction can make 
the electric control of the magnetization dynamics dras- 
tically more efficient [8] and give rise to giant spin mo- 
tive force [9] , greatly raising the technological prospect of 
spintronic devices based on unprecedented strong mutual 
coupling of magnetization dynamics and electric current. 

Compared to other spintronic phenomena such as gi- 
ant magnetoresistance and spin Hall effect [10, 11], theo- 
ries of the Rashba interaction in ferromagnets are largely 
phenomenological [2-4, 7-9], based on free electron- like 
Hamiltonian such as 



T~i — 7^0 L He 



(i) 



where Ho = p 2 /2m gives the dynamics of conduction 
electrons in the parabolic approximation, % ex = — Jcr ■ n 
is the ferromagnetic exchange coupling between the con- 
duction electron spin cr and the localized moment n, and 



Hr = ckr(<t x p) • z is the Rashba spin-orbit coupling 
interaction, with the surface normal to a Rashba ferro- 
magnet taken as the z-axis. In these theories the Rashba 
parameter cvr is a purely phenomenological parameter, 
and it remains unclear what controls the strength of «r. 
A naive estimation of an based on the relativistic ef- 
fective magnetic field in free electron picture and con- 
sequent Zeeman energy results in ~ 10~ 4 eVA, which 
is several orders of magnitude smaller than experimental 
values [6, 12, 13]. In light of the technological prospect 
of Rashba ferromagnets, not to mention the theoreti- 
cal importance of the problem itself, we try to estab- 
lish the Rashba interaction in ferromagnets on a micro- 
scopic level. In recent papers some of the present authors 
showed that the Rashba interaction in non-magnetic 
bands is essentially a multi-orbital phenomenon [14, 15]. 
A similar idea was proposed earlier in Rcf. [16]. As the 
magnetically polarized bands also typically exhibit multi- 
orbital character with most of the d-orbitals involved 
in the band structure, the multi-orbital scheme may be 
brought to bear on the magnetic system as well. 



Construction of a microscopic Hamiltonian is done for 
tig d-orbitals (d xy ,d xy ,d zx ) where magnetism typically 
occurs. We assume a square lattice and introduce the 
tight-binding Hamiltonian H t2g for electron hopping to 
nearest neighbor sites with Slater-Koster parameters t\ 
and ti for a- and tt- hopping of the d-orbitals. The inver- 
sion symmetry breaking (ISB) about the xy-plane gener- 



2 



ates additional hopping terms [17], 

HlSB = l^W^xy^Ci+V,™,* + c lxy,*^+£,yz,°- + h - c A 
i . a 

~ 7 ^ ] i C i.xy.cr C i-y,zx,a + c i^ X y. a c i-i,yz,a ^.C.], (2) 



dictated by the ISB parameter 7. Here i denotes the 
atomic site and Ci tXy / yz / ZXta - is the electron annihilation 
operator of the d X y/dy Z /d zx orbital with the spin z- 
componcnt a. In addition, magnetic exchange H cyi and 
atomic spin-orbit interaction (SOI) Hamiltonian i/soi 
arc introduced as 



Hsoi 



(3) 



where Cj, 



(c f 



y,cr> z,y2,CT' z,za;,(T 



), and L is the 3x3 



matrix consisting of the expectation values of the atomic 
orbital angular momentum (OAM) within the ti g orbital 
space. This completes the microscopic model we shall 
study now: H TB = H t . 2g + ff ISB + H cx + H S oi [18]. 

Numerically determined energy bands of Htb are 
shown in Fig. 1. The result indicates interesting depen- 
dence of the band dispersion on the relative direction of n 
with respect to the Bloch momentum k. When n is cither 
parallel or anti-parallel to k we do not find any changes 
in the dispersion [Fig. 1(b)]. On the other hand, when 
n is orthogonal to k, a clear displacement of the band is 
observed [Fig. 1(a)] in the opposite directions according 
to whether n || +z x k or n || — zxk . The amount of 
the displacement grows linearly with k near the T-point 
(k = 0). We remark that Fig. 1 is obtained in a situation 
where J is larger than the other energy scales so that 
the electron spin direction is essentially parallel or anti- 
parallel to n. Then the behavior of the energy bands in 
Fig. 1 can be summarized as the energy shift caused by 
the Rashba interaction (cr x k) • z in the presence of H cx , 
demonstrating that the phenomcnological Hamiltonian 
T~i [Eq. (1)] can indeed be derived from the microscopic 
one, Htb ■ This dependence of the energy bands on the 
n direction has been utilized in the experiment [19] to 
verify the Rashba interaction at Gd surfaces. 

The n-dependent band shift noted in Fig. 1 van- 
ishes when 7 = 0. To capture the essential role of 
the ISB parameter in establishing the Rashba inter- 
action we focus on iJisB = St -ffiSB(k), written in 
momentum space through Ck,<r = A^ -1 / 2 J^. e lk ' ri Gj )Cr 
[iV=number of atomic sites]. Near the T-point we have 

#iSB(k) = -2701, So- Ck.aClj ' k x z) c Kv, where a L 
is the lattice constant. Thus L tends to align along 
the ±k x z direction. Wc then combine this ten- 
dency with iJsoi = Sk^soi(k), where iJ SO i(k) = 
^s°Y<acr' C tal L ' ( cr )^(r'] c k, -' • Roughly speaking, the 



combined effect of -ffisB (k) and i?soi (k) is to replace L 
in TJsoi(k) by k x I, producing the Rashba interaction 
proportional to (k x z) ■ a. 
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FIG. 1: (color online) (a) Energy bands for the tight-binding 
Hamiltonian Htb with (ti, £2, 7, A BO , J) = (1, 0.2, 0.5, 0.5, 3.0). 
Red solid and green dashed lines are obtained for the unit 
magnetization vector n directed along +x and —x directions, 
respectively. Majority (M) and minority (m) spin bands are 
labeled with superscripts. Signs of the Rashba parameters for 
the bands are indicated. The band a = 3 shows no Rashba 
effect as indicated by the absence of band shift. Energy dis- 
persion along the k y direction are shown, (b) No discernible 
difference in the band structure occurs along the k y direction 
for magnetization n directed along +y and —y directions, (c) 
Chiral OAM pattern for a — 1 majority band with A so = 0. 
Other parameters are the same as in (a). 

This heuristic argument can be verified by a pertur- 
bation calculation. In the absence of -ffsoi the bands 
separate into three majority (cr || n) and three minority 
(cr || — n) cigenstates |k, a) <8> |<x ■ n = ±) (a = 1, 2, 3) of 



H 



-H K 



t 2g i-iJisB + H CX . Within such band basis the spin-orbit 
Hamiltonian i?soi(k) becomes 



3 

HS0l(k) = ^ i'l^n^alp^a 

a=l 



(4) 



where ijjk,a 



gives the spinor consisting of 



c k,a,+r 
^Ck,a, — r 

the majority (|| +n) and the minority (|| — n) a-band, 
and (• • ■ ) gives the inter-band matrix elements which are 
ignored. The intra-band 2x2 matrix elements is obtained 
near the T-point as 

H k ,i = -H K 2 = a R (k x z) • tr, H k , 3 = 0, (5) 

where aR = (A so OL)(7/fi) is the Rashba parameter. 
Near the T-point the eigenstatcs also carry the chiral 
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OAM [15] 



27, 



(k, l|L|k, 1) = -(k, 2|L|k, 2) ~ — k x z, (6) 

h 

whereas the eigenstates in the third band does not have 
OAM [15]. (Results over the full Brillouin zone are shown 
in Fig. 1(c) for a = 1 band.) 

Putting all the pieces together and returning back to 
the real space, we arrive at the effective Hamiltonian 
truncated within each majority-minority spin-pair band 
a as H TB ~ J2 a H a with 



H a = / d 2 r tft 



2m* 



Jcr ■ n+a R>a (tT X p)-z) *r,ai(7) 



where m* is the effective mass of the band a. Because of 
the large value of J /an, the inter-band matrix elements 
in Eq. (4) only makes corrections which are second order 
in k and are ignored in Eq. (7). Now the Rashba pa- 
rameter becomes band-dependent, an,i = ckr = — ckr.,2 
and even aR.,3 = 0! Such non-trivial band dependence of 
the Rashba parameter cannot be understood from phe- 
nomenological consideration alone. As a result when the 
magnetization direction n is reversed from +x to —x, the 
corresponding shift in the band structure should occur in 
the opposite directions along the fc y -axis for a = 1 and 
a = 2 and not at all for a = 3, as is indeed the case; Fig. 
1(a). For n = +y —> —y the shift in the fc^-direction is 
minimal; Fig. 1(b). 

A few remarks are in order. Firstly, upon compar- 
ing Eq. (5) to Eq. (6) we deduce that the sign of 
the Rashba term in each majority-minority pair band 
a is correlated with the chirality of OAM for that band 
pair. The OAM non-carrying band a — 3, in turn, docs 
not have the Rashba interaction. Both nonzero chiral 
OAM and Rashba interaction are consequences of multi- 
orbital character of the band and the inversion asym- 
metry. With the recently developed circular-dichroism 
angle-resolved photoemission technique one can indepen- 
dently probe the orbital chirality of the bands [20] and 
deduce the sign of the Rashba parameter for the spe- 
cific band in question. Secondly, the connection between 
i?iSB m Eq. (2) and the combination L • (k x z) is nat- 
ural in view of symmetry, because the combination vio- 
lates the inversion symmetry along z-axis and can thus 
arise generically when ISB occurs. In other models de- 
scribing, for instance, p-electrons moving in 2D triangu- 
lar (relevant for Bi) or square lattice, it can be shown 
that the additional hopping Hamiltonian [such as -Hisb 
in Eq. (2)] allowed by ISB is also proportional to the 
combination near the T-point. To be more strict, a more 
general Hamiltonian L • [R(k) x z], where R(k) is an 
arbitrary odd function of k consistent with the crystal 
symmetry, also violates the inversion symmetry. In fact, 
it can be verified that H\sb in Eq. (2) is of this form with 
R(k) = — 2j[sm(k x dL)x + sm(k y aT,)y], which reduces to 



— 27<2£k near the T-point. Note that R(k) x I explains 
the OAM pattern in the entire Brillouin zone in Fig. 1 (c) . 
Thirdly, although the derivation leading to Eq. (7) was 
done for small SOI, one can equally well solve the large- 
SOI limit and obtain the same Hamiltonian as Eq. (7) 
in the total spin J = 1/2 basis with aj = 7a/, [15]. 

Band-to-band variations of a R , a imply that experimen- 
tally determined Rashba parameter [6] will be an average 
of an, a in some sense. To clarify the nature of the averag- 
ing, we derive the Rashba-induccd adiabatic STT [3] from 
our approach. To this end we now regard n in Eq. (7) as 
position-dependent and introduce a r-dependent SU(2) 
rotation of the operator ^ r = U r ip r to adjust the spin of 
the conduction electron to localized spin, U}(a ■ n)U r = 
a z [21]. We focus on effects of one particular majority- 
minority band pair (subscript a is omitted) and return to 
the multiple-band case later. Approximating the emerg- 
ing SU(2) gauge potential —ill^d^Ur by their diagonal 
components, —iU$d^U r ~ a^a z [21], one arrives at the 
effective action S, including the Berry phase action for 
spin n and the spin Hamiltonian H n , as (/=spin size) 



S = -2Ial 2 I dtd 2 v a Q - / dtH n 



dtd 2 r ij)l[idt — aoCz]V'r 

2 



d td 2 r $ [P + aa 
2m 



— J dtd 2 v (qr(ii x [pcz + a]) • z) ip T 



(8) 



Landau-Lifshitz equation of motion for n follows straight- 
forwardly, 



(ia L 2 + p s )d t n + n x 



Sn 



+ (f • V)n + 2ma R n x (f x z) = 0. (9) 

The spin density takes on the usual form p s = ip'a z ip, 
while the spin current is modified, 

— (4'^ <T z [p^]~[p^}a z ifj) + ( — +a R z x n) ^ip, 
Ira \m / 



(10) 



due to the Rashba interaction. The last term on the 
l.h.s. of Eq. (9) is the adiabatic Rashba-induced STT 
discussed in Refs. [2-4] 

For multiple bands crossing the Fermi level, the terms 
involving the spin current in the above equation are gen- 
eralized to (restoring band index a) 



(11) 



In light of the experimental condition which can only ac- 
cess the total current, J s , it becomes evident that the 
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experimentally measured qr 
ramctcr aR. ff defined as 



is 



the effective Rashba pa- Rashba-type interaction becomes symmetry-allowed. 



a R,aia = a R,eff J S 



(12) 



Due to the band-to-band fluctuation of aR. a , aR, e ff is & 
sensitive function of the band structure and the position 
of the chemical potential. Rashba-induced non-adiabatic 
STT discovered recently [7, 8] cannot be derived without 
further including impurity terms in the effective action. 

So far our discussion has implicitly assumed magnetic 
band which itself is subject to ISB 7 and atomic spin- 
orbit coupling A so . In conventional 3d fcrromagnets, an 
estimate of the Rashba parameter q;r ~ (A so flL)(7/ii) 
yields at most 0.1 eVA since A so <C lOOmeV typical of 3d 
transition metals. In contrast, experimental values in the 
range of 1-2 eVA has been reported in 3d magnetic thin 
films covered with heavy-element layer [6] , which suggests 
that mixing of atomic orbitals could result in hybridized 
bands with both strong magnetism and large Rashba in- 
teraction [22]. In non-magnetic systems, the enhance- 
ment of ckr through alloying has been reported [23]. 

The idea is readily confirmed by the first-principles 
calculation with, for example, one Co layer and three 
Pt layers forming a heterostructure where the inversion 
symmetry is naturally broken [24]. Isolating the eight 
bands around the Fermi level, band structures are calcu- 
lated with the spin quantization axis constrained along 
+X and —x directions, respectively. The energy disper- 
sion measured along the fcj,-direction for the two cases 
are plotted in Fig. 2(a). The two bands at the top of 
the figure are nearly non-magnetic and experience little 
Rashba shift in response to the reversal of spin axis. For 
the remaining six bands which are magnetic, we deduce 
significant Rashba parameters |q;r = 1 — 2c VA from 
the band shifts, in excellent quantitative agreement with 
recent experimental values [6]. Composition analysis re- 
vealed substantial fraction of both Co and Pt atoms for 
all the bands in support of the existence of new hy- 
bridized bands with both Rashba and magnetic charac- 
ters. Similar to Fig. 1, the Rashba parameter is heavily 
band-specific [Fig. 2(a)]. In comparison a single mag- 
netic Co layer yields almost no band shift as seen in 
Fig. 2(b). Even imposing perpendicular electric field 
of 1 eV/A failed to improve the Rashba effect for the Co 
monolayer. As another comparison we carried out band 
structure calculations with one Co layer sandwiched be- 
tween two Pt layers on either side to ensure inversion 
symmetry. Here we again find almost no Rashba shift of 
the band. It can be concluded, then, that twofold condi- 
tions are to be fulfilled for realizing magnetic bands with 
large Rashba coupling. One is the hybridization of mag- 
netic and non-magnetic heavy atom orbitals to ensure 
the effective mixing of magnetism and spin-orbit inter- 
action. The second is the inversion symmetry breaking 
such as available at surfaces and interfaces to ensure that 
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FIG. 2: (color online) (a) Band structure of lCo-3Pt layers 
with the magnetization forced along +x (red full curves) and 
—x (green dashed curves) directions. Vertical arrows connect 
bands related by Rashba shift. The band associated with or- 
ange arrow is minority band while others with blue arrows are 
majority bands. Upward (downward) arrows imply positive 
(negative) Rashba parameter. The shift direction is reversed 
for the minority band, (b) ICo layer, (c) 2Pt-lCo-2Pt layer. 
There is little Rashba shift observed in (b) and (c) between 
the two magnetization directions. 

Essential aspects of the mixing can be understood 
by analysis of the tight-binding Hamiltonian, H = 
ff™ + H M + H NM ~ M , consisting of heavy non-magnetic 
(NM) bands with Rashba interaction, light magnetic 
(M) band, and their coupling (NM-M), respectively. 
Through the coupling, each Hamiltonian H a in the non- 
magnetic Rashba band acquires a second-order correction 
in the hybridization parameter t a between non-magnetic 
a-band and the magnetic band, 



(13) 



H a ^ y^V'k.q ( £ k,a + «R,q(q" X k) • Z 
k 

\ta? 



(£k,o — Ek.Mj" 

Correction to the magnetic band is 



Jcr ■ n 4 a 



H 



M 



\t a \ C*R,a 
(ek.a — £k,M)' 



;(er X k) ■ 2)lpk,M- ( 14 ) 



Apart from some energy dependencies in the denomina- 
tor, both Hamiltonians (13) and (14) have remarkable 
resemblance to the effective Hamiltonian previously de- 
rived, Eq. (7). Based on insights from these perturbativc 
calculations we claim the bands shown in Fig. 2(a) are 
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at once magnetic and carry a substantial Rashba param- 
eter by virtue of the heavy mixing of magnetic and non- 
magnetic orbitals. The magnetization dynamics of the 
hybridized bands is also governed by Eqs. (7) through 
(12). 

The present investigation showed how magnetic bands 
with substantial Rashba interaction can arise in mixed 
magnetic and heavy-element structures. In contrast to 
standard "relativistic" picture of Rashba effect, multi- 
orbital nature and atomic hybridization are crucial fac- 
tors in forming the enhanced Rashba magnetic band. 
Currently available experimental values of or are very 
much scattered [6, 12, 25] even for apparently similar 
magnetic layer structures. We speculate that the sensi- 
tive band dependence of the effective Rashba parameter 
as discussed in Eq. (12) may be partially responsible for 
this. Away from the T-point the simple linear-k depen- 
dence of the Rashba interaction should be replaced by 
some odd nonlinear function of k, which will modify the 
angular dependence of Rashba-induced STT as well. 
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